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Abstract 

The picosecond carrier dynamics in a closely packed Si-nanodisk (Si-ND) array with ultrathin potential barrier 
fabricated by neutral beam etching using bio-nano-templates was investigated by time-resolved 
photoluminescence (PL). The PL decay curves were analyzed as a function of photon energy by the global fitting 
method. We show three spectral components with different decay times, where the systematic energy differences 
of the spectral peaks are clarified: 2.03 eV for the fastest decaying component with a decay time r = 40 ps, 2.02 eV 
for r= 300 ps, and 2.00 eV for r= 1.6 ns. These energy separations ranging from 10 to 30 meV among the emissive 
states can be attributed to the coupling of wavefunctions of carriers between neighboring NDs. 

Keywords: neutral beam etching, silicon nanostructures, bio-nano-templates, time-resolved photoluminescence, 
carrier dynamics 



Background 

Semiconductor nanostructures have been paid much 
attention because of their potential applications to op- 
tical devices with improved performances or novel func- 
tionalities. Quantum effects in the nanostructures are 
essential to realize those optical applications. The 
quantum effects can also affect the carrier dynamics in 
the nanostructures because carrier transfer or separation 
can be understood by a combined effect of the spatial 
distribution of carrier wavefunction with the energy re- 
laxation of carriers. Therefore, the size, density, and 
interspacing of the nanostructure, which strongly modify 
the carrier wavefunction, need to be precisely controlled 
for the device applications. For this purpose, we have 
recently reported Si-nanodisk (Si-ND) arrays fabricated 
by neutral beam etching using bio-nano-templates [1]. 
This method allowed one to provide closely packed 
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high-density Si-ND arrays. These nanostructures were 
difficult to prepare by means of conventional self- 
assembly methods using co-sputtering and annealing. In 
this characteristic nanostructure, wavefunctions of 
photoexcited carriers can spread over neighboring NDs 
because of the periodic alignment of the NDs with ultra- 
thin potential barriers. This spatial overlapping of the 
wavefunction will result in an ultrafast carrier transfer 
among the NDs, which can be induced by the ultrafast 
energy relaxation of the carriers into NDs with lower 
potentials. It is important to point out that superlattice- 
type Si-NDs based on the ND alignment used in this 
study are efficient for future applications to quantum 
dot solar cells with high efficiencies if photoexcited elec- 
tron-hole pairs are immediately separated after the exci- 
tation. In this report, we investigate the carrier dynamics 
induced by the lateral coupling of the carrier wavefunc- 
tions in a single layer of a two-dimensional Si-ND array 
by means of time-resolved photoluminescence (PL). 

Methods 

The Si-ND array was fabricated by damage-free neutral 
beam etching using a bio-template composed of ferritin 
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supramolecules. Details of the fabrication process are 
described elsewhere [1-3]. The ND thickness, diameter, 
and interspacing distance were designed at 8, 10, and 
2 nm, respectively. Optical excitation was performed by 
second-harmonic fs pulses with a wavelength of 400 nm, 
pulse width of 150 fs, and repetition rate of 76 MHz, 
produced using a mode-locked Ti:sapphire laser. The 
time-resolved PL spectra were observed at various 
temperatures and excitation densities using a synchros - 
can streak camera. The time width of the instrumental 
response curve was less than 15 ps, and the time reso- 
lution of 5 ps was estimated after deconvolution with 
the instrumental response. 

Results and discussion 

The time-resolved PL spectra of the Si-ND array at 
250 K with an excitation power of 50 mW (approxi- 
mately 8.4 mj cm" 2 ) are shown in Figure 1. The PL 
emission band appeared at around 2.0 eV. The PL in a 
visible light region was characteristic in many types of Si 
nanostructures. The PL peak energy systematically shifts 
towards the lower energy side as the time elapsed, where 
each PL spectrum at different time regions is fitted by a 
single Gaussian function. From this fitting, the energy 
shift of the PL peak is quantitatively clarified to be about 
30 meV up to 8.0 ns after the pulsed excitation. No 
significant time-dependent change in the spectral shape 
or width was observed. 

The observed PL cannot be attributed to the indirect 
bandgap emission affected by the quantum confinement 
effect, which was often reported in studies using small 
Si nanocrystals with diameters of 2 to 3 nm [4,5], be- 
cause this observed emission energy which is observed 
in the present NDs is remarkably higher than the value 
of 1.7 eV which is reported in small Si nanocrystals 
with an averaged diameter of 2.5 nm [4]. Moreover, life- 
times of the observed PL, which is in the range of 
40 ps to 1.6 ns and will be discussed later, are much 
shorter than those of the microsecond scale, character- 
istic of the indirect bandgap recombination of carriers. 
A possible origin of the present PL is the emissive sur- 
face states which are weakly located around the inter- 
face of Si/Si0 2 . There are also several reports for 
surface-related emissions in the visible region, and it 
has been confirmed by PL measurements of some sam- 
ples with different surface treatments [6]. It should be 
noted that the relatively strong PL with the very fast 
decay times implies the good optical qualities of these 
Si-NDs without defects. The other presumable inter- 
pretation is the quasi-direct gap emission. Dhara et al. 
have reported the PL emission with the wavelength of 
about 600 nm and decay times of several nanoseconds 
[7]. They assigned this PL to the quasi-direct bandgap 
emission in heavily strained Si-NCs. 
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Figure 1 Time-resolved PL spectra of the Si-ND array at 250 K. 

Solid black lines show the fitting result with a Gaussian function. 



The spectral linewidths of single Si nanocrystals were 
reported to be 100 meV or more [8,9], which were also 
dependent on the fabrication method and surface condi- 
tions. In our case, the size of the Si ND was precisely 
controlled by the diameter of the Fe core which formed 
in a cavity of the ferritin molecule. The size uniformity 
of 8% was calculated from the statistical analysis of SEM 
images [1]. The full width at half maximum of the PL 
band observed in this study is in fact 200 meV, which is 
a not-so-large value as compared with the above results 
for single nanocrystals. These considerations reveal that 
inhomogeneous spectral broadening and the resultant 
spectral diffusion due to the size distribution of Si NDs 
are not significant in the present study. In addition to 
this, the rise component of the PL decay curve was not 



Kiba et al. Nanoscale Research Letters 2012, 7:587 
http://www.nanoscalereslett.eom/content/7/1/587 



Page 3 of 5 



apparently observed even in the lower energy side of the 
PL spectra. The absence of the rise component indi- 
cates that the time-dependent spectral shift observed in 
Figure 1 is not caused by the spectral diffusion, which is 
a time-dependent red shift of the PL energy originating 
from the carrier migration among high-density NDs for 
seeking lower potential minima. 

Figure 2 shows typical PL decay curves at different PL 
energy regions in the Si-ND array, where the time- 
resolved PL spectra were divided into 20 energy regions. 
The time-sliced curve of PL intensity in each energy re- 
gion was fitted using a triple-exponential function taking 
convolution of the instrumental response function into 
account. The temporal evolution of each PL decay could 
be expressed by neither single and double exponential 
functions nor stretched exponential function. At least, 
triple exponential components were required to obtain 
the best fit. It was also difficult to determine unique par- 
ameter sets when we increased the number of compo- 
nent, though we could reproduce the experimental data 
well. Therefore, we applied the global fitting analysis 
with linked decay times in order to avoid ambiguities of 
the fitting parameters. For this analysis, all decay curves 
were fitted simultaneously using three universal decay 
time constants, and we obtained individual amplitudes 
of three decaying components for each energy region. 
Then, we could reproduce the PL component spectra 



associated to the decay times by plotting the time- 
integrated intensities as a function of the photon energy 
for the three decaying components. The result is shown 
in Figure 3, where we have identified the three PL decay- 
ing components with different time constants t x = 1.6 ns, 
r 2 = 300 ps, and r 3 = 40 ps, respectively. 

We shall provide a possible explanation for the exist- 
ence of the three PL components with different decay 
times. We assigned these three decaying components to 
three different types of emissive sites in the Si-ND array 
from the disk-density and excitation-power dependences 
of the PL decay time and intensity [10]. First, the emis- 
sion component with the slowest decay time t x on the 
order of 1 ns was interpreted by electron-hole pairs 
localized at individual NDs because this PL component 
was observed to be dominant when we excited lower- 
density NDs with the disk interspacing larger than 
40 nm. This trend was independent of the excitation 
density, which suggested the absence of the excited state. 
In contrast, the two faster PL components with r 2 and 
r 3 appeared only when we excited the high-density ND 
array with relatively high excitation powers. Therefore, 
these faster emissions were understood by the recom- 
bination of electron-hole pairs which are not strongly 
localized in each ND, where each wavefunction of the 
carrier spreads over neighboring NDs. The fastest PL 
component with r 3 was particularly attributed to the 
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Figure 2 PL time profiles and fitting results, (a) PL time profiles of the Si-ND array at different energy regions. Solid lines are the result of 
global fitting analysis, (b, c, d) Typical fitting results of the PL time profile using a triple exponential function, where the PL time profile is 
deconvoluted with an instrumental response function. A black line shows a fitting calculation, and each decaying component is shown by a 
narrow line. The upper panel shows a residual difference between the experimental and simulated values. 
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Figure 3 Decay-associated PL spectra. The spectra were obtained 
from the global fit analysis with linked time constants of a triple 
exponential function expressing the PL decay curve. Solid lines show 
fitting curves using Gaussian functions. A peak energy for each PL 
component is also indicated by an arrow. 



recombination involving the electron transfer among 
the NDs which are induced by the periodic regular 
alignment of the ND separated by ultrathin potential 
barriers. In other words, the fastest PL was quenched 
by the electron transfer. 

As can be seen in Figure 3, the decay-associated com- 
ponent PL spectra show the different energy peaks. The 
peak energy of the component PL spectrum associated 
to the slowest decaying component {I x with ri) is mean- 
ingfully lower than those of the faster two components 
(I 2 with r 2 , and I 3 with r 3 ). The energy difference be- 
tween Ii and I 2 is 20 meV, and that between l x and I 3 is 
30 meV, which are quantitatively obtained from the fit- 
tings for the decay-associated component PL spectra 
with Gaussian functions. These energy differences mean 
the energy separations among the energy levels of indi- 
vidual emissive electronic states discussed above. The 
lowest energy state associated to the slowest PL compo- 
nent (Ix with Ti) originates from the electron-hole pairs 
that are weakly localized in surface states in individual 
NDs, which corresponds to the ground state in the 
present high-density ND system, and the wavefunctions 
of carriers are localized within one ND. Higher-energy 
levels with faster decay times (I 2 with r 2 , and I 3 with r 3 ) 
are excited states in the ND system, and the coupling of 
carrier wavefunctions among NDs can form these 
excited states. Therefore, we consider that the obtained 
values of 20 to 30 meV of the energy differences among 
the decay-associated PL components can be explained 
by the energy differences between the ground state 
and the excited states in the present coupled Si-ND 
array structure. We have also confirmed this result from 
the temperature dependences of the decay-associated 
PL intensities. The intensities of these three spectral 



components were measured with varying temperatures, 
and the results were fitted by calculations taking activa- 
tion energies between each energy level as well as the 
barrier level into account. The differences of the activa- 
tion energy also indicate the energy differences of the 
lower two states, and the obtained value is 55 meV be- 
tween the slowest component (ti) and the middle com- 
ponent (r 2 ), which agrees with the results obtained from 
the global fitting method. Further study to test the 
energetic structure experimentally elucidated here in 
the Si-ND system is needed using three-dimensional 
calculations of the carrier wavefunctions based on the 
actual ND array structure. 

Conclusions 

The carrier dynamics in the high-density Si-ND array 
fabricated by neutral beam etching using bio-nano- 
templates was studied by time-resolved PL, where the 
PL decay curves were analyzed at various photon ener- 
gies by the global fitting method. The obtained decay- 
associated spectral components show clear differences of 
the peak energy of 20 to 30 meV between the slowest 
component (ti) and the faster two components (r 2 , r 3 ). 
These energy differences for the emissive states can be 
attributed to the different coupling states of the carrier 
wavefunctions in the present high-density ND system. 
The ultrafast carrier dynamics elucidated from the time- 
resolved PL analysis demonstrates that high-density 
semiconductor nanostructures precisely designed and fab- 
ricated by neutral beam etching using bio-nano-templates 
are promising for optical applications based on quantum 
effects such as quantum dot lasing and photovoltaics. 
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